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Abstract 

This thesis deals with instabilities found in the infinite, 
collisionless, homogeneous plasma with no external fields. The plasma 
has random motion and contains an arbitrary number of streams, each 
stream having its own mass, temperature and streaming velocity. The 
region of validity of Landau damping is found to include disturbances 
with phase velocities greater than two and one half times the thermal 
velocity in one-component plasmas. It is concluded that marginal 
inverse Landau damping cannot explain the two stream instability. 
Based on analysis of the dispersion relation and numerical solutions 
on the computer, instability is found to increase when the largest 
streaming velocity is increased or the largest Debye number is de- 


creased in a multi-component plasma. 
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Chapter 1 


INTRODUCTION 
1.1 Purpose and Scope 
Particle wave interactions in a hot plasma have been under 
intense investigation ever since Landau! theoretically predicted that 
spontaneous oscillations are damped even though there are no collisions. 


é that such microscopic collisionless damping is caused 


It is well known 
by particles absorbing energy from the disturbing wave when they are 
travelling slightly slower than its phase velocity. Conversely parti- 
cles may give energy to the wave by travelling slightly faster than the 
phase velocity. This effect tends to cause the wave to grow. Normally 
both these effects are present and it is a matter of which one dominates 
that determines whether the wave damps or grows. 

Such microscopic phenomena occur when there is a distribution in 
random motion. On the other hand macroscopic bunching occurs when drift 
motion is present. This instability is caused by density fluctuations 
in the plasma. (see section 2.8). This thesis will investigate 
plasmas where both random motion and drift motion are present. One of 
the objectives of the thesis will be to determine the conditions under 


which each interaction dominates. As well, a method for determining the 


stability of a multi-component plasma will be given. 
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The model to be considered is a hot multi-component plasma, 
each component with its own velocity distribution, density, tempera- 
ture, mass, and streaming velocity. The plasma is infinite and homo- 
geneous throughout, with no external electric or magnetic fields. 

Only two types of free energy sources to feed the instability 
will be considered. The first, the macroscopic source, is derived 
from the drift kinetic energy associated with each of the beams of 
particles in the plasma. The beams are streaming with velocities 
relative to each other. This is the macroscopic source of free energy. 
The second, the microscopic source, is derived from the energy re- 
leased when the initial velocity distribution relaxes toward thermal 
equilibrium. 

This model does not include bounded magnetically confined lab- 
oratory plasmas, but places special emphasis on the importance of the 
plasma parameters such as Debye length, temperature, and mass, and 
their effect on the stability of the plasma. This approach will no 
doubt assist in the physical interpretation of instabilities as well 
as allowing the results to be extrapolated to more complicated cases. 

Basically there are two mathematical models to describe a plasma, 
the hydrodynamic and kinetic formulations. Since the hydrodynamic 
equations are in reality velocity moments derived from the kinetic 
equations, they do not describe the detailed behaviour of the gas but 
only the macroscopic phenomena occuring in the plasma. The kinetic 
formulations, however, are able to account for the microscopic fluctu- 
ations as well. Landau damping was first discovered! using one of the 


basic kinetic equations, the collisionless Boltzmann equation. This 
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equation together with Maxwell's equations will also be the 
mathematical basis for this thesis. 

The literature involving the combination of microscopic and 
macroscopic phenomena in plasmas, for the type of model described 
above has mainly dealt with the cases in which a distribution in ran- 
dom velocity is assumed in combination with at most two streams. 
Macroscopic phenomena for the same model, have, however been dealt 
with for the multi-stream case. In this regard, James and Vermeulen3 
assumed a distribution in drift velocity in studying macroscopic 
phenomena in a microparticle plasma. This thesis will therefore try 
to bridge this gap in the literature, studying both microscopic and 


macroscopic phenomena using the collisionless Boltzmann equation. 
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Chapter 2 
DERIVATION OF THE DISPERSION RELATION 


Ee basie eCuat1ons 


Various derivations of the dispersion relation for longitudinal 
and transwerse waves in a warm plasma are given in Schmidt”, 
Bernstein” and Tanenbaum®. The present treatment will make a brief 
review of the derivation of the dispersion relation. This will allow 
the introduction of all the important parameters and will show the 
limitations involved. Since the type of plasma to be considered is a 
multi-component gas, each component with its own given density, temper- 
ature, mass and streaming velocity, the Boltzmann equation for each 


component is: 


Vio =[" 0 (21.1) 


where a is the distribution of the o component and is a function of 
position, velocity, and time. 
here dG is the electric charge and m, is the mass of particles of 
the o species and v is the velocity, E the total electric field, B the 
total magnetic field, 
ty 
ot 


is the collision term. 
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Equation (2.1.1) together with Maxwell's equations, given below 


in MKS units, form a self-consistent set of equations. 


(2ele2) 


where p(r,t) is the charge density and Eo the permittivity of free 


Space. 


Also, 
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Where U(7,t) is the current density and c the velocity of light. 


The distribution function will be represented by 
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Where ee is the equilibrium velocity distribution of each species and 
only velocity dependent and tile is the perturbed velocity distribution 
of each species. The assumption that Wap ie will be made, since 

only linear interactions will be considered. Also, in equilibrium it 


will be assumed that space charge neutrality holds. 


Because there are no external electric and magnetic fields, these 


will be treated as first order quantities: 
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Before the solution to the above equations is attempted, the 
collision term appearing in (2.1.1) and the zeroth order distribution 


function will be examined. 


2.2 The Velocity Distribution Function 

Distinction should be made between an equilibrium distribution 
and the most probable velocity distribution of all the particles in 
the plasma. When equilibrium is achieved, the velocity distribution 
should no longer be changed by collisions, hence the collision term 
in eauation 2.1.1 is then equal to zero. The most probable distribu- 
tion function is the one for which the Hamiltonian is a minimum. 

It is quite well known® that for a monatomic gas with a given 


temperature, density and mean velocity, both the equilibrium and the 
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most probable velocity distribution is the Maxwellian: 


2 
a) seu 
fo (vy) = =o 
Oo ie) as (2neely) 
eo} 
where liv 
kel, 
gees m, Zn e) 


a is the thermal velocity, kp Boltzmann's constant and ge the 
temperature, n is the density while aaete U, is the drift velocity. 
In this thesis U will be assumed to be in the x direction. For the 
plasma model being considered, the Maxwellian distribution (equation 
2.2.1), for each species is the most probable distribution function 
Since the Hamiltonian is a minimum? The plasma is not necessarily 
at equilibrium since this is dependent on the collision term. For 
example, if the Boltzmann collision integral’ is assumed, which is 
used for two body collisions, the range of interparticle force laws 
being short, the most probable distribution function (equation 2.2.1) 
is not the equilibrium one since the collision term does not vanish 
for long times. On the other hand, if the Krook collision oe is 


assumed: 
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Where v-~ is a heuristic relaxation time, the collision term vanishes, 


for times longer than yl, 


The Krook model (also called the B.G.K. due 
to Bhatnagar-Gross and Krook) has some serious physical limitations.°® 
For example, every average macroscopic value approaches equilibrium at 
the same rate v7!, Actually, in the collisions between particles of 
two different species the change in the lighter species' momentum is 
more pronounced than the change in the lighter species' kinetic energy, 
hence the equilibriation time for the mean energy of the lighter species 
is much longer. Therefore this collision model is only valid in the 
limit when v approaches zero for a multi-component plasma. The most 
probable distribution function as given by equation (2.2.1) is the 
equilibrium distribution function in the limit of a collisionless 


plasma. The analysis in this thesis will be restricted to the limit 


where the plasma approaches the collisionless state. 


2.3 The Dispersion Relation 


Returningeto equatiOnse2 lsc Uole2ul. 2 and SupSUItutiInGg, 2c. C0 
2.2.3 into them, gives after performing differentiations with respect 


to time and space: 
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Solving equation (2.3.1) for i eegives: 


ee = = 5 
“Egl8g Fs, x i 


es z (a8) 
m. pole a G vor o xy. k) 
Solving equation (2.3.4) for B, yields: 
eS 
B= x Fy (en) 


Equations (2.3.6) and (2.3.7) are now substituted in equation 


(2.3.5) to obtain the dispersion relation: 
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The transverse and longitudinal oscillations are uncoupled if there 

is no external magnetic field, or inhomogeneouties in density or tem- 

perature. | Therefore longitudinal waves (KI |E,) and transverse waves 

(KIE,) will be considered separately. The perbutation will be taken 

in the x direction. Hence k=ky. By noting that the left hand side 
2,0 


of equation (2.3.8) has only a transverse component: -ic‘k Ee 


5) 
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and by dividing through by iw, the following relation can be obtained 


from equation (2.3.8): 
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debye length and the components of c, are defined by 
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2.4 Longitudinal Waves 


The x-component of equation (2.3.9) yields: 
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where 1 = une Vector in x-direction. 


To evaluate the integral in equation (2.4.1), the quantity 
ctu, is substituted for V. (see equation (2.2.1)). By performing the 


dot product as well, the following equation is obtained: 
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with j = y and z are identically equal to zero because the integrand is 


odd with respect to the y and z co-ordinates. 


Hence there are only two non-zero terms in equation (250 2) 
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By dividing through by EY the following equation is obtained: 
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Notice that only the streaming in the x direction contributes to 
longitudinal waves, and that the second and third integrals are each 
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equal to 7 a. 
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In the limit of a collisionless plasma 
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If both the numerator and the denominator of the integrand in equation 
(2.4.6) are multiplied by ain z) and if it is noted that the integral 


over all odd functions of p vanishes, the following is obtained: 
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The integral can be simplified in the following manner. 


G(s) is introduced with: 


(2047) 
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By multiplying equation (2.4.10) by e°“o and integrating from s=0 to 


S=1, the following is obtained: 
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Now from equation (2.4.9) G(0) = —< ea (a4 13) 


Since z is complex, the poles p = + Zz, are located anywhere in the 
complex plane. For every Z., one pole p is below and one is above the 


real axis. Therefore the theory of residues may be used to show that 
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G(0) = 92 i 


In arriving at this result, the first part of the path of integration 
is chosen along the real axis, and the second part is chosen along the 
contour formed by a semi-circle about the upper half plane. 


Hence, solving equation (2.4.12) for G(1): 
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Gli) isea function of Ze and is renamed Z(2 ). Hence 
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By going back to equation (2.4.8) the dispersion relation for longitu- 


dinal waves may now be put in terms of Z{2z) and so 
az 22) (2.4.16) 


2.5 Transverse Waves 
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Since EY = 0 for transverse waves and Uy is in the x direction, 


equation (2.5.1)may be written for the y component: 
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A similar equation for the z component exists. 
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In arriving at these expressions (equations (2.5.2)and (2.5.3) )an 


analysis similar to that in Section 2.4 is used. 


By multiplying the integrand in equation (2.5.3) by Sean) ez) 


this equation becomes: 
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It will be convenient to replace w/k by Von Since V_, is complex, it 


ph 
will be called ‘complex phase velocity’ in this thesis. 
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Hence equation (2.5.5) becomes: 
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2.6 Qualitative discussion of the Dispersion Relations 
The solution of equations (2.4.16) and (2.5.6.) is difficult with- 
out some knowledge of the behavious of zAca)). This function has been 


tabulated by Fried and Conte. / By using the relationship’ 
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-2(1 + Ze Z(z_)), C2621) 
equation (2.4.16) may be expressed as follows: 
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Only instabilities where k is real are examined (universal instabilities). 
The method of solution of equations (2.6.2) and (2.5.6) (in this thesis 
referred to as ‘method A'), as outlined by Tanenbaum? is straightfor- 


ward: 


1. A suitable complex phase velocity w/k is chosen. 

2s Z,» one for each component is found from equation (2.4.5). 

3. k is found from equation (2.6.2) for longitudinal waves and 
from equation (2.5.6) for transverse waves. 

4, Since k is now known, w can be found from the value of the 
phase velocity assumed in step I. | 

5. By trial and error, values for w/k can be obtained such that k 
is real. If the imaginary part of w is found to be positive, 


the solution is unstable, while a negative imaginary implies 


a stable solution. 
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From step 5 above, it follows that all complex phase velocities which 
result in unstable solutions lie in the upper half Voh plane. Since 
the sign of the imaginary part of Voh is the same as the sign of the 
imaginary part of Z,> the poles Zz, of the integrals in (2.4.6) and (2.5.4) 
also will lie in the upper half Zz plane when unstable solutions 
exist. Conversely stable solutions correspond to values of Voh and Z 5 
in the lower half of the Voh and Ze planes respectively. 

Since the pole Zo corresponds to the oth stream in the plasma, 
there will be as many poles Zz, as there are streams for each Von The 
location of each pole Z in the Zo plane is a function of the Von of the 
disturbance, and of the streaming and thermal velocity of the oth stream 
as can be seen from equation (2.4.5). For every case considered in this 


_ thesis, a reference frame will be chosen such that the stream with the 


lowest drift velocity (called the first stream o=1) has u,=0. 


Consequently 
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SOmulal Z4 will always lie in the same quadrant in the complex plane as 
the Voh of the disturbance. 
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Every pole z, will contribute a real part, kn (Re Lia e)) or 
Ka (Re DL a) to the dispersion relations and an imaginary part 
kA (Im TAZ) or ke (Im Liz) ) For k to be real, «2 must be positive 
and real. 
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Stream on the stability of the plasma. 


The real and imaginary parts of the functions A\z7) and Ne) 
have at most 3 maxima or minima and are always finite for constant 


imaginary zy in the upper half plane, as can be seen from the tables 


if 


given by Fried and Conte’. As a increases the humps are much less 


pronounced. On the other hand, when Ze is in the lower half plane, 


the functions oscillate rapidly for |z. 


if 
Tee > Za . The peak 


amplitudes approach infinity as |z. 


incr 
areal creases for constant z., 


eal" 

In the lower half plane because of the many zeroes of Im Ne 
and Im ON there are over 50 zero's for constant imag(z.)=-9), 
there will always be values of the complex phase velocity in the lower 
half Voh plane for which Ko is real and positive. Consequently, if no 
values of Von can be found in the upper half Yoh plane for which K2 is 
real and positive, it may be concluded the plasma is stable. 


For transverse waves, without solving equation (2.5.6) explicity, 


jt will now be shown that this model predicts that no instabilities 


exist for transverse waves, assuming oA a With this approximation, 


equation (2.5.6) becomes, using equation (2.6.1): 
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When solutions exist, Re (k7)> 0 and for unstable solutions, Im(V,,)>0. 


However, the right hand side of equation (2.6.4) is always negative 


SincewkerZ.\zZ.)2-2 for Im(V,,)70 from the tables of Fried and Conte’. 
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: is real and positive 


Since no values of Von! << c exist for which k 
when Von is in the upper half Voh plane, this model predicts that trans- 
verse disturbances in plasmas are stable when Baral ser Under these 
restrictions, instabilities for transverse waves are possible, however, 
when velocity distributions are anisotropic® or when there is an 

applied magnetic field’, These cases will not be studied in this thesis. 


The solution for longitudinal waves (equation 2.6.2) will now be 


considered in more detail. 
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2./ Effect of Parameters on Distribution Function 

The upper half plane is given in Figure 2. 

Region N is where Re{ Z'(z,)| is negative 

Region P is where Re|Z'(z.)| is positive 

he right-hand plane is where Im[Z"(2,)] is negative 

The left-hand plane is where In[ 2" (2) | is positive 

For unstable solutions to exist, there must exist at least two 
poles, not all of them in the same quadrant in the upper-half 7A plane. 


Otherwise Z ee ImZ'(Z.,) # 0 and then alee # 0 as can be seen from 


Oo 

equation (2.6.2). At least one pole must lie in region P so that Ke 
can be positive. Poles in region N will be referred to as stabilizing 
and poles in region P as unstabilizing. 


By using equation (2.4.5) the poles Z, may be expressed in terms 


of V h as follows: 
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n a, 
To test for instability, given all parameters (eo ays u,) of 


each stream in the plasma, every possible point Voh in the upper Voh 
plane must be checked until solutions are found by method A in Section 
2.6. As can be seen from Fried and Conte’, Re Li Ze) is even and 
ImZ'(z.) isvodd (form constant Im(z,). A reference frame such that u,=0 
and Uns Uz ++ > 0 is chosen. Then for any ay chosen in the upper 


left-hand V_, plane, the Zz, 5 also lie in the upper left-hand Z5 plane 


ph 
and so no instabilities exists. Hence Voh need only be considered in 


the upper right-hand Voh plane. 
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The choice of Von can be restricted even further by the following 
method. The largest magnitude in streaming velocity will be denoted 
by Une Then only the Re(V4) eoue need be considered, since for 
Re(V 4) Zod all poles Zo lie in the right-hand Z5 plane and no un- 
stable solutions can exist, as explained earlier. According to the 
tables of Fried and Conte’, tf Im(z.) eeaIRE\ 2 i. the pole z. will lie 
in region N. Hence from equation (2.7.1), if Im(V 4) rolls all poles 2 


will lie in region N. In summary, only V_, such that Re(V..) < u. and 


ph ph m 


Im(V,) <u, need be checked in the upper right-hand plane for unstable 


solutions. 


It will now be shown, that, for constant Im(V iy) there will be 


at least one Re(V.,) for which Im (k>) = 0. A reference frame is chosen 


such that u,=0 and U5, Ug... 2 0 as above. If there is at least one 
Stream # such that u. > 0, when Re(V..) = 0, 4 k “Im[Z" (z )] >0 as can 
71 ph Bete O 
: 5 Z : 
be seen from Figure 2. Moreover, when Re (Vy) Sle x kG Im| Z (z,)]<0. 


Hence, for any given Im(V 4)» there will be at least one Uy such that 


) If and 


te 5 . 
y kOe Im|Z (z.)] = Qsand this Uh will be denoted by (V e: 


pi 2 
only if 2 kno Re| Z (z.)] & nen 


ph 


Shy , will there be a solution. 


ph oh) s 
Next it will be shown that, if there exists a solution, there wil] 
; Z 
be at least one th (=v) for which ? ko Re| Z oy = 0. When there 


| 
exists a solution, © ae Re Z'(2,) 2 UETOr. Im(V.) <u, When Im(V 1) 


? | 
Sales all ae will lie in region N and so 2 kg Re| Z (z.)] < 02 Hence, 


if the plasma is unstable at any O<Im(V,)<u, there will be at least 


. 2 ee 
one (Von), = Vo for which y oe Re|Z (z,)|=0. 
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Consider any plasma where all the parameters u 


been chosen or determined except Us 


chosen high enough, the plasma will be unstable. 


Now, u. is chosen such that: 
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Ve: Re(z.) a 


i (all o) 


where Re(z.) = 


a. 
Re(z.) > Re(z 


The subscript s refers to that component for which 
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It will be shown that if Un as 


(aefe2) 


From the graphs of Fried and Conte’, there exists a Zp such that 


Re|Z'(2,)| is positive for values of [Re(z_)| > |Re(z,,) | for constant 


Im(z,) albco, 


At Uap S72 heehee and Re(z,)>0, 
a Re(z.) + u. - u 
os a S Som 2 

Re(z,) = a Re(z.) 

Im(z.) = Im(z.) 

Re(z_) a oRe\ 7 teu ase Goins ai 

uv a a 

m m 


The contribution to (2.6.2) from Bs 


Re|Z'(z,)| and In| Z(z,)] are small. 


(25753) 


(2.7.4) 


(287 25)) 


is very small since both 
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Also, 
~ eae Im(z 
Im(z__) = "s a (2.7.6) 
an 
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Hence all poles lie in region P, so that Retk eo. and all poles 
other than Zz, }ie in the right hand plane so that Im(k“)<0. 


Moreover at en =e at Re(z,) + a, Im(z,) and Re(z,)<0 


tee) (2.7.9) 
0,8 Ons 


Since Re(z_ _) is very large, it gives a negligible contribution. 


Rea ne a ) (2.7.10) 
an Un 7 Un 

Re(z_ ) = ae 2454 + a = Re(z,) (Zeal) 

Im(z,_) = Im(z,) (Jee) 


2)50, Since all poles 


a) 


Again all poles lie in region P so that Re(k 


except Z, give a negligible contribution, Im(k™)>0. 


Since for one value of Yoh? mers 120) and for another value 


2 


2 
Im(k~)>0, there exists at least one value of Von such that Im(k~)=0 


while Re(u +z, a,,)>Re(V,)>Retu.+Z, a.) since u, is so much higher. 
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For any of these values of Von? Re(k-)>0 Since all the poles stayed in 
region P. This proves an instability exists if Un is large enough. 


Uren will be defined as the minimum bes for which instabilities exist. 


Ueeeen Sn i > 
( eves ot necessarily >> Zea) 


Any change in the 2 j 
y g parameters aos Kno ee We etc. which cause an 
increase in Ua will be referred to as stabilizing since the range 
over which u, can vary for an unstable solution has been reduced. 


Likewise any change in these parameters which cause a decrease in eas 


Wilisbewreterced sto. as UnStabilizing. 
The definitions of important and negligible poles will be 
explained next. In a plasma where the number of components is greater 


than two, some poles Zo might give negligible contributions compared 


to others. For instance, this may happen when Iz, 1>>|z,| (kj 00 kp 52) 


in which case the contribution of Zz, may be considered negligible. A 


: : . 2 
pole Zz. 1s called important when any change in ass ky; Or u. causes a 


marked increase or decrease in u The terms important and negligible 


min* 
poles will be used in conjunction with the term stabilizing pole which 


may be used to describe a pole Z. when, Re[Z'(z;)]| <0. As explained be- 


fore, when ReZ'(z.)<0, Z is in region N. 


2 re 
How the parameters Ueda Kno and Te affect the position of the 


poles and the stability of the plasma will now be discussed. If Un is 


increased, the range of Im(V oy) and Re(V 4) to be tested will be 


: Z 
increased, resulting in a greater range for Im(V 4) for which Re(k~)>0. 
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Usually this is an unstabilizing change. 


tO 


a> Ao > [Re(z.) | < 0.924 


and consequently Z. will be stabilizing. On the other hand, if 
a; << Up» Iz | will be quite large for Monn 4, | ~ O(u,) and therefore 
makes a negligible contribution. For the small range where 


\V -u, | ~ O(a.) the pole makes a significant contribution. 


ph 


The Debye number squared (ky25) is an effective weighting factor. 


Ine the Timrt that dee approaches zero, the pole z; can only make a 


small contribution to Ke, For the case in which oy is increased 


such that 


2 + 
Pimikee 27.) te (tas 
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no solution exists. 


For the case when the pole Zz lies in region N a solution is 


2 
impossible provided ar is large enough to ensure that Zk Re ze} 


<a ince is increased with a large enough increase in ky 2 


min o-, suéh 


an increase of ke is a stabilizing change. Hence an increase in the 


largest Ke can usually be classified as a stabilizing change, 
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Mass, ms is inversely proportional to a, and is independent of 


and a are proportional to k),é and independent of 


2 : 
Koo . Density, le 


ay: Consequently, the deductions made when a, is varied also hold when 


2 2 : 
Ym is varied and similarly, the deductions made when k,“ is varied 
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also hold when i and a are varied. 

Temperature, on the other hand, is a function of both the 
thermal velocity and Debye number. The pole Zz, is inversely 
proportional to ies and ae is inversely proportional to le. If 


ie is increased from the initial temperature 103 to a critical value 


2 2 
lobes such that KDoe <<kp 4 and such that a 0(0,0), the pole Zc will 


give a stabilizing contribution, since it is in the region near 

the origin in Figure 2. But theweighting factor ae, will make 

the contribution from Zee ely small compared to the initial contribution. 
Hence, ie can be increased such that ra makes a negligible contribution. 


It will now be shown that for T. sufficiently small (T_.) the 


os 
contribution to Ke will always increase. ive is assumed sufficiently 


small so that Z os will be sufficiently large that the approximation 


ee =) = - = can be made. This is the asymptotic expansion for 


os 
. 7 ane 2 
Z ee) given by Fried and Conte’. The contribution of Zo to k 


is given by: 
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The contribution to k” from the pole Ze will then be, in terms of 


the initial contribution: 


2 

Pee 
D 

: a (2.7.18) 
bz i 


Hence, for ee sufficiently small, which implies b small, the 
contribution from equation (2.7.18) will be greater than that from 
equation (2.7.13). Therefore, a large decrease in ie will make the 
pole Z, more important. 

Section 2.7 will now be summarized. 

1. Increasing Un, increases instability. 


2. Increasing the largest eae increases stability. 


U 
m . esyip etna 
3. When a. > 5 90g > it causes pole Z. to be stabilizing. 


4, When a 55 Us it causes pole Z. to make a negligible 
contribution. 

5. Increasing T. by a large amount causes z. to make a 
negligible contribution. 

6. Decreasing T, by a large amount causes Z. to make an 
important contribution. 


The summary above describes the usual cases. The various physical 
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mechanisms responsible for these instabilities will be reviewed in 
section 2.8. Then, in Chapter 3, using section 2.7 and the analysis 
of some special cases with the aid of the computer, the regions 
dominated by the different physical mechanisms (bunching and Landau 


damping) responsible for most of the instabilities will be established. 


2.8 Physical Interpretation Review 

Balancing the energy lost by particles to that gained by the 
disturbing waves, and using the dispersion relations of various 
plasma models, the complex interactions can be understood in the 
following way?, 

In a plasma composed of at least two streams, a mechanism 
called bunching occurs. The electric field of a random disturbance 
in one stream causes an acceleration of particles in the second stream 
and a fluctuation in their density. The particles can carry this dis- 
turbance back to the original starting point of the first stream if 
they have streaming motion. If this latter disturbance is in phase 
with the original one, the wave will grow, so that an instability 
occurs. Thus, assuming one component at rest, for every different 
streaming velocity u,, there exists a mode for which the disturbance 
is in phase with the original one. Since the original disturbance is 
random, the different modes corresponding to different beams must be 
added to give the total effect. If the nodes of one mode fall on 
the troughs of the other, destructive interference occurs, called phase 
mixing, lessening the bunching effect. If there are many streams, 


this phase mixing tends to dampen the original disturbance. 
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Another interaction between particles and waves already mentioned 
in section 1.1, involves particles travelling at the phase velocity of 
the disturbing wave in a warm plasma. If there are more particles 
travelling infinitesimally slower than the phase velocity than there 
are travelling infinitesimally faster, the wave will lose more energy 
than it will gain, and damping occurs. If the reverse is true, that 
is, the slope of the velocity distribution is positive at the phase 
velocity, the wave will absorb more energy from the interacting parti- 
cles than it will lose and instability ensues. This type of damping is 
generally known as Landau damping while the instability is sometimes 
called inverse Landau damping. 

Classically“ bunching is considered to occur in cold plasmas, 
and Landau damping in hot plasmas. In Chapter 3 it will be shown 


where the regions of validity for each instability actually are. 
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Chapter 3 


QUANTITATIVE DISCUSSION 


3.1 General 


In order to illustrate the various interactions described in 
2.8, consider the following plasmas. A warm, one component plasma 
with velocity distribution given by the Maxwellian (see equation 
2.2.1) will be called plasma A. Assuming Landau damping dominates, 
all waves will be damped since the slope of the velocity distribution 
is negative everywhere. Section 3.2 will consider this case in 


more detail. 


The cold plasma approximation to plasma A can be described as 
follows: an infinite number of cold beams each with its own velocity, 
but each beam infinitessimally spread apart in velocity space such 
that the distribution of velocities is Maxwellian. If bunching is 
assumed to dominate, this plasma will be stable, since the infinite 
number of streaming velocities will suffer phase mixing, as explained 
IMeseCulony 2.0. 

As a third example consider a plasma composed of two hot streams, 
each stream like plasma A. If inverse Landau damping dominates, the 


plasma will be unstable because part of the slope of the velocity 


distribution is positive. However, section 3.3 will show that even 
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though the slope of the distribution function is positive for a range 
of values of the phase velocity, the dispersion equation (22622 ewe 
not allow growing waves for this range. The physical explanation of 
Landau damping emphasizes only the few particles which are in resonance 
with the disturbing wave, while bunching tends to emphasize the 
collective process of many more particles in the plasma. The warm 
dubble-humped distribution function will serve as an example to 


illustrate this point in more detail in section 3.3. 


As a last example, consider two cold streams far apart in 
velocity space. This plasma will experience bunching and be unstable 


as described in section 2.8. 


A thorough investigation of the dispersion relation (equation 
2.6.2) using the distribution function (2.2.1) will help to determine 


the regions of validity of bunching and inverse Landau damping. 


3.2 The "One-Component" Plasma 


For a one component plasma, ions forming a fixed background, 
there exists only one pole. From Fig. 2, for only one pole Zy> 
Re[Z'(z,) ]<@ when Im[Z'(z,)]=0. Hence only in the lower half plane 
of dy will there be any poles for which Im[Z'(z,)] is zero, and 
Need positive as explained in Chapter 2. This is a stable 
configuration and all waves will be damped when the component has a 
Maxwellian velocity distribution. The problem has been solved and 


discussed in various ways and at various lengths in textbooks and 
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papers For the least damped solution and for phase velocities 


greater than 2% to 3 times the thermal velocity, the damping factor 


is very close to 
aeewe 
L =i eK a 
= 1 u) 4 & P 
pater se ne) te eee oe Ane o ( 
Ws Bae ly) 
i 3 3° 
Where w., the imaginary part of w, is the damping factor and w,, the 


seh p 
plasma frequency. In terms of the distribution function, given in 


equation (2.2.1), and after integration over the transverse velocity 


components, the damping factor (3.2.1) is proportional to 


W) °F (ua /k) 


2 ; (3272) 


nk 


That is, the damping factor is proportional to the slope of 
the distribution function in this case. Therefore, for one component 
plasmas Landau damping prevails for disturbances with high phase 
velocities. As can be seen from equation (3.2.1), for phase velocities 
greater than approximately three times the thermal velocities, the 
actual damping factor becomes very small. Because there are very few 
particles interacting with the wave, these plasma oscillations can be 
considered as free modess Hence, Landau damping is considered an 
exponentially small effect. 

However, for phase velocities smaller than and of the order of 


the thermal velocity, the damping factor no longer is proportional to 
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the slope of the distribution function. In these regions, the damping 
factor is very large’. The region of validity of Landau damping, 
the interaction between the particles and the wave at the phase 


velocity, includes, therefore, those disturbances with V_.> 2.5a,. 


ph ] 

The large damping factor for Von < 2.5a, may be caused by the phase 
mixing of the disturbances associated with bunching. However, this 
can only be determined when an analysis of phase mixing for warm 
plasmas has been carried out, and at present, to the writer's know- 
ledge such a study has not been made. Qualitatively any Maxwellian 
velocity distribution should exhibit phase mixing. 

Because bunching is caused by density fluctuations, it depends 


on Debye length. Moreover, because of its dissipative effects, high 


thermal velocities will diminish bunching. 


3.3 Two Stream Instability 
The following study of the two stream instability will show 
the regions of validity of inverse Landau damping and bunching. These 
regions will be established by studying the distribution relation 
jr Ok 


will be investigated. If the two stream instability is in the 


(equation 2.2.1). First, the marginal stability case Ore 


regime of marginal inverse Landau damping, equal numbers of particles 


are travelling slightly slower and slightly faster than the Von velocity. 


Hence, 
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Where ene is the one dimensional distribution function of the o 
component of the plasma. If solutions can be found for which 


yi er # 0, these solutions will not be in the regime of marginal 
oO 


inverse Landau damping. 


By substituting equation (2.4.15) into equation (2.4.16) the 


following dispersion relation is obtained 


2 aoe z 
ae 2 aez _ 6 
(6) 


Since the poles are located close to the real axis, the real part 
of equation (3.3.2) is approximately given by the first two terms on the 
right hand side, and the imaginary part is given by the last term on the 
right hand side. The imaginary term can be expressed in terms of the 


distribution function (equation 2.2.1) and so, 


z 
(0) 
Paes. ie meee Z 
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Ronee Bin gee lay e2zh J. g diem ieee (Gac-3) 


where the transverse velocities have been integrated out of the three 
; ; 2 ; 
dimensional distribution function. Since Im(k™) = 0 for solutions 


to exist, as explained in Method A of Section 2.6 and 2.7 
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Using equation (2.3.10) this condition becomes: 


al 
Bee NouaDn! = (3.3.4) 
O Mm, 


From equation (3.3.4) it can be seen that there will be no solu- 
tions which lie in the regime of marginal inverse Landau damping 
(equation 3.3.1) except for equal charge squared, mass ratio of both 
Streams. This is why in general an explanation different from Landau 
damping should be sought for the warm two stream marginal instability. 
In this regime bunching would appear to dominate. 

The above might seem to imply that the regime for inverse Landau 
damping includes two streams with equal Debye length. However, by 
analyzing this case for complex Von? jt will be shown that this case 
too, belongs to the regime where bunching is probably valid, rather 
than inverse Landau damping. 

The cases which correspond to the regime of inverse Landau 


damping are those for which 


> ft [Re(Vop) | >0 (erets) 
©) 


By using the dispersion relation (2.6.2) it will be shown that, 
even when 
: (3.05.6) 
nf [Re an) | 0 


instabilities may still exist. 
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For a two component plasma with equal Debye numbers and thermal 


velocities, the dispersion relation becomes: 


ye (z,) i zlze) (Cor oa7) 


where ky is the Debye number. The poles zy and Z, are complex. The 
Streaming velocity Uy 1s set equal to zero, as done in section 2.7. 


From equation (2.4.5): 
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Where Wy is the real part of the frequency, and subscripts have been 
dropped from a and u. W; must be positive for an unstable solution to 
exist. When Im(z_) is fixed, the symmetry properties of La Nen) are 
such that 1m[Z' (z,) }is an odd function of Re(Z,,) and Re{ Z'(z,)] is an even 
function. Therefore, for the imaginary part of the right-hand side of 


equation (3.3.7) to disappear 
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By solving for Re(V,,) = =. equation (3.3.9) becomes: 
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Using equation (2.2.1) and (3.3.10), the sum of the derivative of 


the distribution functions at the phase velocity gives 


2 
S| = ate 


ae 
gC sa Elect) 


Equation (3.3.6) holds rather than equation (3.3.5) showing that this 
case is not in the regime of inverse Landau damping, but indicating it 
is in the regime of bunching. 

Since it has been established that the two equal stream insta- 
bility does not lie in the regime of inverse Landau damping, it may be 
assumed that most multicomponent plasmas are likewise not dominated by 
inverse Landau damping if they have instabilities. 

It may also be assumed that disturbances with (Von - u.) rs 3a, 
will be damped because they are in the regime of classical Landau 
damping. Because of the above, efforts will be concentrated on para- 
meters such as thermal and streaming velocity and Debye number. 

The functional relationships between the streaming velocity, 
plasma frequency, and thermal velocity for the marginal two stream 
instability will now be discussed in more detail. Let the stream with 
the lowest ky be denoted by i, the stream with the highest Ky DV ih 
For each stream there exists one pole 2Z,, (k=i,j) where Zz, is given by 
egudeion (2.7.4). solving tor Voh in terms of Zz. and substituting in 


the equation for Zs the following expression is obtained: 


S : 
sSWijr, 2 16 
(LEE A) 0 ; : . | 


atid dad ontwoie (2.2.4) norteune nat? hdrey, ebTad (d.€ Pan 


fH pakieotbha! suc, paTqnsh, webnss sevsvii to swipsy ott oF Ton as a 
. ourdoagd To sith sag ate at 

<Bianiinaavde {apo awt Sad Tanith bere’ cares nase. ped 7! ete 
ad an. de aeotonsh Yeodhd s27QVNt. 70-501 aS 544 ar oF! ton ly 
vi berantingh ton gerwai |! ava Renae ty spe nectieu ed fae jeom MmAd fs . 
sersetVaaiant aved qth) }i aad vabrsd Si 

a Ale dae iW agoibiweyalt Sod Vetues os) pate) i oa 


Wabi Teorerata te suieoreah) Wi ans Yam Seuho o! senna. ty 


iy 


wide wel ’eteniaesnde. sd | UW ertdtys oveds err avuguad < 
wStli\in, aviea0 Dre Yiroolay whines se file rages £8, e 
ITIOFSY balingntt2 oi _Bswial et enor IES tenotesnyt 
maond 2 ow! lohieienosit or yireotsy tanya ors pEmsoooNT Me 
ASWinee ae int so.) Chtadsh aun nf bacevonit @eh wom Pw vi 
Nd gt testy rHisds Ari sige) as »t ya nena ad a 
NaonevEp 2t 8 so5itw CL ote} us ater 410 ateiwe sissy / cord 
mi ‘pilfus (tae bile ys vi ing" ng anivtor abt. ie 
+ p re 4 aa akan; xen ni mu 
ww : 


43 
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Since Im[Z' (2) is odd for constant Im(z_), ReZ. )and Re(Z;) must 
have opposite signs for Im(k2) to equal zero in equation (2.6.2). 
Since Im(a;2;)=ImV=Im(a,z,), the imaginary terms in (3.3.12) cancel. 


Therefore equation (3.3.12) can be written as follows: 


ae Z bit ed 2) y= lu, - uU (gooe hs) 


where Zz. = Re|z,| and z, = Re|z, | 


The quantities a, u are determined by each stream while Z. and 
Zz, are restricted to such quantities that Im(k*) = 0 and Re(k°) eeOnain 
equation (2.6.2). Figure 3 shows z. as a function of zy) for various 
values of ky a/kp;- This graph can be used to determine whether warm 
two stream plasmas are unstable or not. If Ze and Z) exist such that 
equation (3.3.13) holds, the plasma is unstable. The minimum relative 
difference streaming velocity required for instability can be deter- 
mined by choosing the quantities Zoo 2 such that the slope of the Z. 


vs. Z, curve is equal to -a,/a, from d fu, =e = 0s 


| 
ds 
This section will be concluded by studying a multi-stream plasma 


which can be reduced to an equivalent two stream plasma. 

In a research note Gautene et al report various machine calcu- 
lations from the Saha equations to obtain realistic densities of the 
ionic species for given electron densities and temperature. The para- 


meters given in that note will be used as an example for solving the 
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dispersion relation (equation 2.6.2). The function ae) given in 
equation (2.4.15) was programmed on the I1.B.M. 360 computer (see 
Appendix ). By using method A as outlined in section 2.6, the onset 
of instability was determined by taking Zo real for various values of 
the streaming velocity for each component. This again was programmed 
on the 1.B.M. 360 computer. The various results obtained for the in- 
Stability of this plasma are summarized in the following tables. 

T refers to temperature 

N refers to the number of ion and electron streams 

U refers to the case when all the ions streams are at 'rest' 
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ion 
Ul refers to the case when only the first ion stream o=1 is at 
rest and all other ion streams are at U.: Similarly U2 refers 


to the case when only the second ion stream o =2 is at rest and 


all other ion streams are at UL: 


The explanation of U3, U4, U5 and U6 are all similar to U1 and U2. 
U, refers to the minimum streaming velocity of the election stream 


required for instability in any of the above mentioned cases. 
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TABLE 3.3.2 Instability of a Multi-Component Plasma 
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TABLE 3.3.3 Instability of a Multi-Component Plasma 
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TABLE 3.3.4 Instability of a Multi-Component Plasma 
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TABLE 3.3.5 Instability of a Multi-Component Plasma 
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TABLE 3.3.7 Instability of a Multi-Component Plasma 
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TABLE 3.3.8 Instability af a Multi-Component Plasma 
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VABLES 3.0.10 Instability of a Multi-Component Plasma 
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TABLE 3.3.1) Instability of a Multi-Component Plasma 
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TABLE 3.3.13 Instability of a Multi-Component Plasma 
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TABLE 3.3.14 Instability of a Multi-Component Plasma 
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These results will be examined by reducing these cases to the 
two stream case. 

When all the ion streams are at rest (case U), the dispersion 
relationship (equation 2.6.2) is the same as if there were only two 
streams, with thermal velocities ay and ay respectively and Debye 


number kot and kyo respectively. In terms of the quantities from case 


U, these new thermal velocities and Debye numbers are: 


ape ks 


ion 
Ae es a, 
Ke 2 i ke 
o=] 2 
kyo = Kpe 


Figure 4 shows the marginal instability for this case in terms 
of ky i/kpe vs U,/a,- 

When only one ion stream is at rest (cases U1-U6) the dispersion 
relation can no longer be used to reduce this to the two stream case 
in the same way as was done with case U. 

The thermal velocities of the ion streams at the same streaming 
velocity as the electron stream are not equal to the electron stream's 


thermal velocity. However, since the thermal velocities of the ion 


streams are much less than the minimum velocity required for instability 


(U,)> the ion streams which are not at rest may have a neglible effect 
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as was discussed in section 2.7. If we ignore those ion streams, the 
cases Ul-U6 may again be reduced to the two stream case. In terms of 


the quantities from the cases U1-U6, 


Ea a wus 


Ko1 = kno where o refers only to that ion stream which is at 


rest. 


Figure 5 shows the marginal instability in terms of ko1/kpe Usa oar 
for the U1-U6 cases for which U,/a, lies in the same range as in case U 
(fig.4). Since the ion streams were not ignored when Ue was calculated 
and both these graphs (fig. 4 and 5) show approximately the same curves 
for marginal instability it would appear justified to ignore the ion 
streams with the same streaming velocity as the electron stream. 

Figure 6 shows two curves for the marginal instability in terms 


of kp 1/kpo NES U,/aq: The ratio a. /a, for the cases U, U1-U6 was 


on 
1/172. The lower curve summarizes the U1-U6 cases for which U/a, lie 
between 1.0 and 4.0. The upper curve is shown for comparison. It shows 
the marginal instability for the two stream case with different Debye 
numbers and streaming velocities, but equal thermal velocities. Table 
3.3.16 show the values used to plot the figure. U; refers to the 
minimum streaming velocity required for instability of the stream with 


the largest Debye number. Stability is to the left and under the curve 


for each case. 


} , 


sy smoayte not seodd, prane! ae 


wane 
Yo cever fl «=. f2ro MaSTte owl ae Sasuste 


* 
, ot 
} 
er BP detnw abate r aGAy 6 1) 21ST 
a)” oo - ' 
ah a 
i 
a 
' 
: Fe a 
my ii BY 7 iia \" 1 tla a bigs gi! a ie yom oy! 2 2 slp - 
. j i ' 4 _ 
: a 7 , @ " 
w 2ert. el UW dod Aot ebeoe GUI Bid We 
, , . 
i \ : | | : Vig nor ii ; 
‘ és 
* | i 
LAL } PHN» hn 4 (¢ sad oy *) inne to S2o/1 
it a y tid nr df 
ive # i- ; 4} 2 f at Ai] { aA yin } 72 Sta? art) 
wed F me wet wer 
‘ j .1 i i ' wi f Vek 3 


PY a Nge Wat Mie TT ea . zysh ania eveuy vewal off — SVEM 

( \.& bre 0.1 saat 

ew ain , rr oT Bon Oy Tv hese ! tL yt! af ad . BTE ai og NM a}4} 
» _ 

quh0 farettih diiw e253 maerte ows off vot yiPbiaetant fentqvem sag 


sjdet s@sitipofay laniniy [rigs tt ,e2*) soley grfemeitz bos 8Gee 
| ; * 

* ra : . a ; 

oft ad erste) (aught) sid Salq: 8 hyde amithgy ort wore Gk se 

7 7 

Atiw moat? ofa tu Va f rt, factoal pas + dt m4 Foote panbnanss, 


a? 


ays 5il} "Stn BnA< “ 45] Hit? Hj 2F Xz! hag ny 3 vain Gi aap | 
7 a 


ie : 7 “oes a. As * 
Ay 


7 
! 


64 


AyL_iqgezs jeurbuey g aunbl4 
Soyag 


SJAUND JHL YJONN OANV 1331 3H1 OL S! ALlldavis 


rota 2 


3VAUD SHT RAIGUY. GMA. 7331 BHP OT 2! ¥Tistsate 


~ 


gir frduse hala. 2 empty 


65 


OV 


AyLLLqeys Leurbury 9 aunbly 
OF Or 


Ihde satiecdsaNNe GN Vela sie Oles| 


Ol 


ALIN@VLS 


(8)Zdy dy 


3yVaU> 3AT g20vu CHA TIAd SAT OT el YTINGATe 


a 
ays 


eaviidess (intpre 3 -sver4 


UABLERS +3) t6GTistabilitysapea Two-Component Plasma 


A5x105 


Sinese: 


0.54 


Sie 


8 


4.00 


oo) 10s 


3 


66 


| me 
aa ad 
: ' ‘ oe 7 " 
6ua6i ansnomod-owl 6 10 | i thee g 


~ 


207.0 o,f 


67 


3.4 Multi-component plasmas 


This section will explore the problem of a multi-component 
plasma with one stream having a high thermal velocity, using computer 
calculations. 

As was shown in Section 2.7, a pole Zz. of stream i is stabilizing 
when its thermal velocity is high enough. Further from Section 2.7, 
this effect is increased when the Debye number Kp; is large. In some 
plasmas, therefore, an increase in thermal velocity causes an other- 
wise unstable plasma to become stable. 


As the first example consider an unstable plasma with the follow- 


ing parameters: 


TABLE 3.4.1 Four-Stream Plasmas 


ai/ay 10 AO Sy aMive 1.0 


U./a. 0 (ease ACEO OTS 28) 


A marginal instability is found at Von’ 21 = 0.843, with poles 


and their contribution to Me as follows. 
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TABLE 3.4.2 Four-Stream Plasmas 


2 2 
j Z. Im(k ); Re(k ); 
1,4 0.843 -0.65 -0.143 
2 -1.4) 0.65 +0.530 
3 -0.0185 0.0 -0.387 


Although the pole with the high thermal velocity (z5) gives a 
negative contribution to Re(k*), the Debye number k53 is not large 
enough relative to Kno to stabilize the plasma. The contribution 
of pole Z3 to I mk?) is almost zero. The functions of poles zy and 
Z4 is to cancel the rather large contribution z, gives to Im(k’). 
The function of pole Z5 is to cancel the negative contribution Z3 
and Z 4 gives to Re(k°). 


The second example shows a case where a plasma is stable. 


TABLE 3.4.3 Four-Stream Plasmas 


1 ] 2 3 ad 
LANNY 0.141 1.0 0.25] 519 
a, /a, 1.00 1.00 1.00 172 
U /a;(max)0.0 4.0 0.0 0.023 


As was discussed in Section 2.7, only for phase velocities 


< 4.0 can there by any instability. For this range, the 
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pole with the high thermal velocity, Za: has a value of approximately 
-0.02 and so the contribution of this pole.to ie is almost zero, 
and to Re(k*) iS approximately -0.54. The maximum contribution at 
any IV nn / ay [<4-0 to Re(k°) of all the poles is denoted by Re(k°), 

and is shown below with the appropriate Von! 34° 


TABLE 3.4.4 Four-Stream Plasmas 


ie Re(k*) Von! 34 
1 1.48  +0.011 1.48 

2 21.48 «40.55 2.52 

3 1.48 ~—-+0.0284 1.48 

A -0.02 -0.54 0.0-4.0 


There is a small range, around Von! 33 ~ 2.5 when Re(k-)>0. The 
contribution of the poles to Tm ko) is denoted by Im(k*), and shown 


below with Von! 43 Cy aoe 


TABLE 3.4.5 Four-Stream Plasmas 


2 
i Z. Im(k F 
] tspaye - 0.00 
2 -1].48 0555 
3 2257 -0.00 
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: 2 
Since Im(k") # 0.0 when Re(k“)>0, the plasma is stable. The computer 


calculations confirmed there was no Von for which this plasma is un- 


stable. 


In the last two examples, the stream with the highest thermal 


velocity had a small streaming velocity relative to its own thermal 
Pa 
a. 

; 
stream has a larger relative streaming velocity. 


velocity =Q . The next example shows an unstable plasma where this 


TABLE 3.4.6 Three-Stream Plasmas 
Streams: 1 2 3 
ky i/Kpo 0.998 1:0 0.914 
a./ay 1.0 1.0 i72 
U-/a; O7416 0 0.505 


A marginal instability is found at Von/aq7t- 71. The poles and 


their contribution to Re(k¢) and Im(k?) are given in the following table. 


TABLE 3.4.7 Three-Stream Plasmas 


: 2 Z 

j Z. Im(k ) Re(k ); 
1 #2 20 -0.85 0.547 
2 ye! -0.30 0.52 

3 -0.50 bas -0.875 


Two streams are needed to balance the contribution to Im(k2) 
made by the stream with the high thermal velocity. This is in contrast 


to the previous examples where this contribution is almost zero. 
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Another difference with the previous examples is the effect on the 
stability of the plasma if the stream with the high thermal velocity is 
removed. In the previous examples, this removal results in an unstable 
plasma with a different Von’ 24 for marginal instability. However, in 
the last example, the removal results in a stable plasma. 

The next example shows a plasma where the value of the highest 


a./a, is less than 3.0 compared to 172 in the previous examples. 


TABLE 3.4.8 Three-Stream Plasmas 


Streams: 1 Ze = 

kp 4/Kpo Oe fant hale 0.944 
ai/ay mae) 1.0 2206 
u./a, 0.0 See 0.0 


A marginal instability is found at Von’ 44 = 2.17. The poles and 


their contributions to Re(k°) and im(k°) are given in the following 


table. 


TABLE 3.4.9 Three-Stream Plasmas 


: 2 re 

j Z. Im(k ). Re(k ), 
i MAY -0.03 O2175 
2 -1.14 1.04 Ose 

5 P05 = 1502 Oe15 


It is worthwhile to compare Table 3.4.8 with Table 3.4.1. The 


values for the streams with a./ay = 1.0 are similar while the 
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value of A/a, in Table 3.4.1 is much higher than in Table 3.4.8. The 
effect of this difference is found by comparing Table 3.4.2 with 


Table 3.4.9. It is seen that the pole z, at marginal instability is 


: 
almost zero in Table 3.4.2 and large in Table 3.4.9. The result is 
that this latter pole does not make a negative contribution to 

Re(k*). 
With the next example, the value of A/a) is a little higher than in 


Tabilews 24 eo. 


TABLE 3.4.10 Three-Stream Plasmas 


Streams ] 2 3 
kya /Kpe O2741 10 0.944 
a /ay 1.0 ese) Dgeye 
u,/a; 0.0 Sea 0.0 


A marginal instability is found at Von! 2 = 2.41. The poles 


and their contributions to Reta) and Im(k>) are given in Table 3.4.11. 


TABLE 3.4.11 Three-Stream Plasmas 


: 2 (a 

71 Z. Im(k ); Re(k ); 
] Bab -0.02 O214 

Z -0.90 1,42 Cait 

$) 0.825 -1.40 -0.25 


Every value in Table 3.4.10 is the same as in Table 3.4.8 except 


the ratio A/a, As can be seen in Table 3.4.11, the contribution of 


An instability is found even though the Debye number is large. 
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Paes 
Z2 to Re(k ) is negative. The pcles zy and 25 contribute positively 
2 : 
to Re(k"), cancelling the stabilizing contribution of pole Z3. 
The next and last example shows a stable plasma. All parameters 


are again the same as in Table 3.4.8 except the value a/a,. 


TABLE 3.4.12 Three-Stream Plasmas 
Streams 1 2 3 
kp i/Kpo 0.741 1.0 0.944 

a5/a, 1.0 1.0 2.98 

u./a; 0.0 S2aL™ SOC0 


The computer calculations showed the olasma to be stable. To 
illustrate this, the following two approximations are made to obtain 


a rough estimate of z, when Im(k2) is zero. 


= kyo 


2. The contribution of pole Z, to Im(k2) is negligible. 
By using equation (2.4.5), the above approximations, and the 
values of the parameters from Table 3.4.12, the following expressions 


are obtained. 


Zo ie (3.4.1) 

Pak 
1 

= - 3.4.2 

Zo = 24 aca ( ) 

aE zi! (3.4.3) 
2.98 

Since Im Z'(z) is even, Zo = -23 for a marginal instability to 


exist. Hence ra 2.48. The poles and their contributions to 
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Im(k”) and Re(k*) are shown in Table 3.4.13 when z, = 2.48. The 


] 
values of kp j that are used in the calculations are the same as in 


hapleros 4et2s 
TABLE 3.4.13 Three-Stream Plasmas 


2 Z 
j z. Im(k™) Re(k™); 
] 2.48 -0.02 +0.19 
2 -0.83 +1.49 -0.30 
3 Oates) -1.33 -0.27 
Re(k°) is no longer positive, hence no instability exists 


and the plasma is stable. 

The last three examples will now be compared. Tables 3.4.8, 
3.4.10 and 3.4.12 are identical except for the value of a3/ay. The 
following table gives a3/a) together with pole zy and Re(k?) 


calculated when Tm(k2) is zero. 


TABLE 3.4.14 Comparisons 


2 
zy a3/a, Re(k™) 
Ext Cea 2.06 +0.64 
Exad BN) 2.92 0.0 
Exes 2.48 2.98 -0.38 
From Table 3.4.14 it is seen that, for these examples Re(k°) is 


decreased when a3/ a) is increased. This confirms that plasmas can be 


stabilized when the value of a. is increased. Section 2.7 showed that 
5 


U 
: : m 
Zz. gives a negative contribution to Re(k~) when a> 7-904 ° When 
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divided by ay this expression becomes: 


jes) 


el: m 
ay ° 0,924 a (3.4.4) 


Using the parameters of the previous three examples, with i=3, and 
U 


_m = 3.31, the expression becomes a2/a, > 3.58 when Z3 gives a negative 
a 
contribution to oe Tables 3.479.. 3: 4.1) and 3.4, 16 show tiat, in 
these cases, Z3 is stabilizing for values of Az/a4 less than that in- 
dicated by equation (3.4.4). The value lies somewhere between 2.06 and 
2.92 for these examples. To actually stabilize this plasma, the value 


A/a, had to be increased to 2.98 which is still less than the value 


in equation (3.4.4) for these examples. 


This concludes the study of particular examples for the multi- 


stream case. 
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Chapter 4 


SUMMARY 


Instabilities in the infinite, collisionless, homogeneous 
plasma with no external fields and with an arbitrary number of streams, 
each stream with its own mass, temperature and streaming velocity 
has been discussed in this thesis. The physical interpretation of 
Landau damping and inverse Landau damping has been investigated and 
it was found that the regimes for which these phenomena were valid 
were not always the classical ones. It was found that, for one 
component plasmas, the region of validity of Landau damping includes 
those disturbances with phase velocity greater than two and one half 
times the thermal velocity. The regime of marginal inverse Landau 
damping was found not to include the two stream case. 

Bunching, when applied to hot plasmas as well as cold, was 
shown not to contradict the plasma dispersion relation. (equation 
2.6.2). Rather, the parameters involved, debye length, thermal and 
streaming velocity, were found to support bunching in hot plasmas. 

By using the dispersion relation an extensive description was 
given how temperature, debye length, streaming and thermal velocity 
changes affect the stability. It was shown that changing any one of 
these parameters does not necessarily mean increasing or decreasing the 
stability always the same way. Rather the value of the other parameters 


of other streams have a direct bearing on the whole plasma stability. 
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Some graphs have been included to aid in determining the 
Stability of a warm two stream plasma, showing the relationships 
between various parameters. A method to solve the dispersion 
relation was reviewed in section 2.6, This method was used to 


program the computer to solve some cases. 


Examples were found where a multi-stream plasma can be made 
equivalent to a two stream plasma. Various examples were discussed 
of multi-component plasmas, with one stream always having a higher 
thermal velocity. Results were compared and found to be in 


agreement with theoretical limits developed in section 2.7. 


The appendix summarizes a method to program the plasma 
dispersion relation in an economical way, combining several known 


mathematical formulations of the dispersion relation. 


77 


by 


we siikatemgrt wt ule oe we 


gg triendiitatos ‘edd es ean aii 2) OW! inian a Weer 
(to hireaiae 1h ont uvfce ot hotgeant * _etahalli hey. 3 
Oj See elw Honsan 24). .Ock Sottoae eerily ie 
2bacn wilok oy ioe oe var, i 


abeniad Mes BuesTS One 29) TiN 6 stile neh ase an 
boviuioriti wie! eolsirde SuwaaY §=.amaaly mosh OWE: 6 OF 
veer er ynived dyowlp wagits She TW OSnEB A Seoqnorgenet=# 
nt st of ettoy) Ob baieqnion 41S aiTuasr “tisely’ 
VS holracr we beqgal send’ ed imi £394 sanpsiit datw: 


cipal AY ively 82. vulleteR saa tarae Pee ae 
nwond fot@ves eC Midhdg vei Tasrinainie ne a gat: atte gota 


nateaion notevaqcth alg +e gnoleelaieie? ror 


FOOTNOTES 


LO 


eS 


Landau, L., Physik Zeits. Sowjetunion 10, 154 (1936). 

Stix, T.H., "The Theory of Plasma Waves", McGraw-Hill Book 
Company, New York, (1962). 

James, C.R. and Vermeulen, F. "A Microparticle Plasma", 
Canadian Journal of Physics, 46, 855(1968). 

Schmidt, G. "Physics of High Temperature Plasmas", Academic 
Press, New York, (1966). 

Bernstein, I.B. "Plasma Oscillations:II Kinetic Theory of Waves 
in Plasmas". Nuclear Fusion, 4, 61(1964). 

Tanenbaum, "Plasma Physics", McGraw-Hill, New York, (1967). 

Fried, 6.D., Conte, S.D., “The Plasma Dispersion Function , 
Academic Press, New York, (1961). 

Buneman, 0. "Dissipation of Currents in Ionized Media", 
Physical Review 155, 506(1959). 

Fried, B.D. and Gould, R.W., "Longitudinal Ion Oscillations in 
a Hot Plasma", The Phsyics of Fluids, 4, 139(1961). 

Gautem, M.S. and Ananandrum, M.N., "Computation of the Composition 
of a Highly Ionized Plasma in Thermodynamic Equilibrium", 
Canadian Journal of Physics, 49, 492(1971). 


Franklin, R.N., Research Notes, Plasma Physics, 10, 806(1968). 


78 


7 


: > 
dneer) tel’ OL aeravivagnoed nae aba? od eed 


= 


Inat Thtteneviae . ‘eevsw sheets Fo Wngent Sat” “wees 


(Se8h) 410 way yngamg) 7 


pe ra 
Bit S'S oot aagoy im AY | sretusyay. his ee ¢2OMsL 
f AR} y c } R F i] | ti? T eter neToane. 
: v i a ube a 
simebn! zanest) Sadan . to tata © 2 ibm? 
jie beh : 
; aoe!) , toy wot tear : 
i.” - Pa 
auel? Yc ait sttontd Ulsanperell rae af§" “2.7 tatetaqys® a 


-P 
- 


ma sory 4 | , “Seteat nt . i 


. 
=) 


e i] 7 
(cS ) MY wav . Vat sweectc vam in Ae 4 | a ee 


ny yi’ nT ; : - 5 «a8 bet 


~ 


* , ji ; r) 
TQS hut wah oaagtd epee ae ; ob 


a ~~ 

| i es! set fy ' ~ 4 4 Wy) re) lia ae ‘nen ; & 
PACT POT ,ceP waiver potowte! i 

Oy; eee Tul) i i yi [ ee cite We ‘ hue e Ad ‘ a.8 bain? ae 


: 
| 1 pes! hn of Th | b 2c2tiyven4 SY a anal 4 TOH 6 
rapsieoqia sag To orter7yeng )” yh airs NARHGAA OWS tum a 
»“quivd! Tid otninaybonist at sues! bestrol vine 6 Fo” 
f 


~ 
((XOPISOW 20 .2oTentt Yo Taawel Wefbaed 7 


(BOCL)A08 OF .2nfeyt? amesl? .2s3o0V" dovesesd y pie - 


ay: 


BIBLIOGRAPHY 


Bernstein, I.B. "Plasma Oscillations: II Kinetic Theory of 
Waves in Plasmas". Nuclear Fusion, 4, 61(1964). 

Buneman, 0. "Dissipation of Currents in Ionized Media", 
Physical Review 115, 506(1959). 

Franklin, R.N., Research Notes, Plasma Physics 10, 806(1968). 

Fried, B.D. and Goud, R-Wo,. “Longitudinal lon Oscillations an 
a Hot Plasma", The Physics of Fluids, 4, 139(1961). 

Frved.3B De. Conte, S.D.. “tne Plasma Dispersion -unctions. 
Academic Press, New York, (1961). 

Gautem, M.S. and Ananandrum, M.N., "Computation of the Composi- 
tion of a Highly Ionized Plasma in Thermodynamic Equilibri- 
um", Canadian Journal of Physics, 49, 492(1971). 

James, C.R. and Vermeulen, F. "A Mircroparticle Plasma", 
Canadian Journal of Physics, 46, 855(1968). 

Landau, L., Physik Zeits. Sowjetunion 10, 154(1936). 

Schmidt, G. "Physics of High Temperature Plasma", Academic Press, 
New York, (1966). 

Stix, T.H., "The Theory of Plasma Waves", McGraw-Hill Book 
Company, New York, (1962). 


Tanenbaum, "Plasma Physics", McGraw-Hill, New York, (1967). 


79 


pe uri nt 
(eet a ate ~ 
sank 403 OL-zatevay cara ane 
af anotdaltto20 nol tanthusioael” a 8 ! ayer 
(eaOt}OL! .& p2htwla Io: evpegn? BAT. y, opi 6 
imran? pokersda rl pent? oar” Me pd 4B bate? a : 
ERE) aia 22ers ghinabeah oh 
-Fangtin) alt 7a of de hitartll chal ibe stitaclenbta 2M. smegued ‘an 


-Hyuit dips cihanybeinad? at smeaTh hecinec teat 6 es 7 


~~ 


(Thel jag ink P geptevel Surtaaual wa ye a a 

-“‘anaal? efajtnsqagjaith yi’ at, Tema Beis: A Py " i 

(RaNL S888 [IM Jeaheus® (to Fenuph qtelaned _- 7 

(RSet her OF KoTMitaiwnd Hes Abeus4. ..d ar 

zen? mabeal |, "sant Sta ragoa! iipeh Yo stayin" a. totodae 


AaORE Te ial 
foot 0) HewbaQoM Masvaw: OneEET. 40 qT aT" ssh wht2 
hPL) aval : 


ANDER) « b¥o't wall CP itewaniiol. ae ronda 


APPENDIX 
A, Plasma Dispersion Relation 
Various expansions exist for the plasma dispersion relation 


Fedudeion (2e4el>)) and 1ts derivative’?! 


Ideally suited for 
machine calculations are those expansions which make use of continued 


fractions. A continued fraction is of the form: 


(A.1) 


Due to convergence problems that are encountered, different 
expansions are required for different regions of the complex plane. 
Petezes xt oy CAgzZ) 


aan are pe Se ele NV 


Z 
Vo Bes 
Z(z)=Frointe* (A.3) 
Where o=2 for y<0 
o = 1 for y =0 (A.4) 
go =.0. for y > 0 


and F is as given by equation (A.1) where 


A, = Z 
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x(3/2-x2 + F) 


Where F in equations A.14 and A.15 are given by equation A.1 where 


Be = - n(n + 4) Aral Ae (Re16) 
B= 2n + 3/2 - x2 je Le (A.17) 


The number of terms needed for convergence in equation A.1 varies 
greatly from about 3 in equation (A.13) to approximately 36 in equation 
(A.7). Since there is some overlap in the four regions mentioned above, 
a choice sometimes exists. The places where overlap exists, expansions 
were used which needed fewer terms for convergence. Expansion A.13 
requires the least number of terms, followed by A.10, A.2, and finally 
A.7. This priority was used to program the plasma dispersion function. 
Values obtained using this program were in agreement with values 


tabulated by Fried and Conte’. 
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